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Abstract 


Time-dependent  finite-difference  solutions  of 
the  two-dimensional  Navier-Stokes  equations,  fully 
coupled  with  the  appropriate  finite-rate  chemical 
kinetic  equations,  are  obtained  for  HF  chemical 
laser  flows.  These  solutions  demonstrate  the  prac- 
ticility  of  Navier-Stokes  solutions  for  chemical 
laser  flows.  Results  are  given  for  steady  flows 
where  large  pressure  gradients  are  calculated  as  a 
natural  part  of  the  Navier-Stokes  solutions.  In 
addition,  the  effect  of  unsteady  fluctuations  in- 
tentionally introduced  at  the  cavity  inlet  are 
studied.  Such  unsteady  mixing  increases  the  laser 
gain  by  more  than  a factor  of  two. 
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Nomenc  1 lature 

speed  of  1 igh* 

mass  fraction  of  species  i 

specific  heat  for  constant  pressure 
specific  heat  for  constant  volume 
binary  diffusion  coefficient  (Djj-Ujj) 

multicomponent  diffusion  coefficient  for 
species  i 

total  internal  energy  of  the  mixture 
internal  energy  of  the  mixture 
internal  energy  of  the  species  k 

vibrational -rotational  interaction 
parameter  - I 

total  enthalpy  of  the  mixture 
width  of  the  nozzle 
enthalpy  of  species  k 

Planck’s  constant 

rotational  quantum  number 

thermal  conductivity  of  the  mixture; 

Boltzmann  constant  in  Eqs.  (11),  (21)- 

(23) 

molecular  weight 
Avogadro's  number 

pressure 

specific  gas  constant 
universal  gas  constant 

vibrational  contribution  to  the  electric 
dipole  moment 
tempera  ture 

independent  time  coordinate 
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column  matrices  based  on  flow  proper- 
ties 

x component  of  the  velocity  vector 
reference  velocity 

y component  of  the  velocity  vector 
independent  spatial  coordinates 
concentration  of  species  i 
small  signal  gain  coefficient 

characteristic  energy  of  interaction 
between  a pair  of  molecules 
mixture  viscosity 

x component  of  diffusion  velocity  for 
species  k 

y component  of  diffusion  velocity  for 
species  k 

collision  diameter  of  molecules,  A 
density  of  the  mixture 
density  of  species  k 

species  production  term  due  to  chemical 
reactions 


I . Introduction 

The  past  decade  has  initiated  the  age  of  hiqh 
energy  lasers,  first  startinq  in  1%6  with  gas- 
dynamic  laser,  and  closely  followed  by  breakthrouqhs 
in  large  chemical  and  electric  discharqe  lasers. 

The  gasdynamic  laser,  which  generates  its  laser 
medium  by  means  of  a vibrational  nonequilibrium  . 
nozzle  expansion,  is  the  subject  of  a recent  book. 
The  electric  discharqe  laser,  which  generates  its 
laser  via  electron-atom  and/or  molecule  collisions 
in  a flow  discharqe,  has  been  reviewed  by  Re  ill . 
The  supersonic  diffusion  chemical  laser,  which 
obtains  a laser  medium  from  the  products  of  chemi- 
cal reaction,  is  nicely  described  in  a review  L 
Warren-’.  All  of  the  above  lasers  involve  the  hiqh 
speed  flow  of  large  amounts  of  gas,  hence,  they  all 
involve  the  realms  of  aerodynamics  and  qasdynamics. 
For  example,  the  HF  or  DF  supersonic  diffusion 
chemical  laser  involves  the  supersonic  mixing  of 
two  dissimilar  streams,  as  shown  in  Fiq.  1. 

Concurrently,  the  discipline  of  computational 
fluid  dynamics  has  become  a third-dimension  in 
aerodynamics,  complementing  both  laboratory  experi- 
ments and  pure  analysis1^"®.  Work  is  advancing  on 
both  numerical  methods  and  applications  to  practi- 
cal engineering  problems  The  present  paper  is  in 
the  latter  vein.  Specifically  it  deals  with  the 
direct  application  of  computational  fluid  dynamics 
to  the  solution  of  chemical  laser  flows. 

The  present  paper  represents  a new,  third 
generation  of  supersonic  diffusion  chemical  laser 
analyses.  First  generation  studies  are  exemplified 
by  the  RESALE  computer  proqram7  (which  assumes  one- 
-dimensional  premixed  flow)  and  the  approximate 
flame-sheet  model  of  Hofland  and  Mirels8*9. 
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Second  generation  studies  involve  more  detailed 
fluid  dynamic  calculations,  such  as  the  boundary- 
-layer  solutions  of  King  and  Mirels'"  and  Tripodi 
et  al".  The  present  investigation  is  a third 
generation  study  which  incorporates  the  solution 
of  the  complete  Navier-Stokes  equations'^  for 
chemical  laser  flows  for  both  the  steady  and  un- 
steady inlet  conditions.  The  full  Navier-Stokes 
equations  are  able  to  model  and  solvethe  compli- 
cated chemically  reacting,  recirculating  and 
separated  flow  regions  at  the  base  of  chemical- 
- laser  nozzles  --  an  important  aspect  that  affects 
chemical-laser  performance,  as  emphasized  by 
Grohs'?.  The  advantage  of  invoking  the  Navier 
Stokes  equations  is  that  such  complicated  separated 
flow  fields  as  well  as  any  lateral  or  longitudinal 
pressure  gradients  induced  by  the  chemical  heat 
release  are  modeled  exactly.  The  apparent  disad- 
vantage is  that  numerical  solutions  of  Navier 
Stokes  solutions  take  a long  computer  time  but  this 
can  be  reduced  considerably  per  run  if  while  doing 
parametric  studies  the  runs  are  mad(  back  to  back 
and  the  final  steady  state  condition  of  one  run  is 
used  as  the  initial  condition  for  the  other  run 
with  slightly  different  inlet  conditions. 

The  purpose  of  the  present  paper  is  to  access 
the  viability  of  Navier  Stokes  solutions  for  both 
the  steady  and  unsteady  chemical  laser  flows,  and  to 
underscuie  advantages  as  well  as  the  present  day 
restrictions  of  such  numerical  solutions.  A 
second  purpose  of  this  paper  is  to  compare  the 
results  obtained  before  and  after  the  kinetics  are 
switched  on,  i.e.  to  compare  cold  and  hot  flows. 

This  numerical  experiment  graphically  demonstrates 
the  effects  of  chemical  reaction  on  the  laser  fluid 
dynamics.  Such  conside  ‘ions  are  important  when 
extrapolating  non-reacti  | supersonic  mixing 
results  obtained  in  the  laboratory  to  the  case  of 
real  chemical  lasers. 

Thirdly  and  more  importantly  it  is  the  purpose 
of  the  present  paper  to  illustrate  the  enhancement 
of  mixing  and  qain  obtained  due  to  intentional 
fluctuations  of  the  inlet  conditions.  To  date, 
the  use  of  the  time-dependent  method  has  been  to 
obtain  steady  chemical  laser  flowfields  as  the 
asymptotic  result  of  large  times;  the  final  steady- 
-state  flow  has  been  the  desired  result,  and  the 
time-dependent  technique  has  been  simply  a means 
to  that  end'4. 34.  However,  the  work  described  in 
this  paper  takes  advantage  of  the  time-dependent 
Navier-Stokes  solutions  to  study  the  transient 
effects  of  unsteady  fluctuations  in  the  flow  pro- 
perties at  the  cavity  inlet  on  overall  chemical 
laser  performance  downstream.  In  particular,  the 
results  show  that  intentional  fluctuations  in  the 
inlet  velocity  vector  can  result  in  increased 
mixing  and  at  least  a factor  of  two  increase  in 
HF  chemical  laser  qain.  This  paper  presents  the 
first  results  of  such  unsteady  mixinq  phenomena 
in  chemical  lasers. 


II.  Physical  Problem 

Consider  the  parallel  supersonic  mixinq  of  a 
stream  of  partially  dissociated  flourine  with  a 
stream  of  diatomic  hydroqen,  both  diluted  to  some 
extend  with  He,  as  shown  in  Fig.  1.  In  the  mixing 
region  downstream  of  the  nozzle  exits,  the  follow- 
ing hypergolic  chemical  reactions  take  place. 


Cold  reaction: 

F ♦ H?  - HF*  ( v)  » H (1) 

Hot  reaction: 

F?  4 H • HF*  (v)  + F (2) 

which  result  in  the  direct  formation  of  vibration- 
ally  excited  HF  as  a reaction  product.  Over  cer- 
tain reqions  of  the  flow,  a total  population  inver- 
sion may  exist,  i.e.  N„F(u+1  }>NHF(v) . where 

N||F(v)  's  the  population  of  the  \ th  vibrational 
level  in  HF.  It  is  a population  inversion  that 
makes  a chemical  laser  work  (see  Refs.  1 and  20 
for  background  on  laser  properties).  In  a chemical 
laser,  this  inversion  may  be  total,  as  previously 
described,  or  partial,  where  the  inversion  is 
carried  by  the  rotational  distributions.' 

In  the  present  investigation,  the  chemically 
reacting  flowfield  in  Fig.  1 is  solved  by  means  of 
the  complete  two-dimentional  Navier-Stokes  equa- 
tions. The  mixinq  is  assumed  to  be  completely 
laminar.  Mu) t icomponent  diffusion  is  treated  in  an 
exact  fashion.  Both  the  hot  and  cold  HF  reactions 
are  included,  and  vibrational  populations  up  to  the 
v - 8 level  in  HF  are  calculated 

The  parallel  mixinq  of  Fiq.  1 was  chosen  for 
solution  because  a)  it  consitutes  a reasonably 
straightforward  test  of  the  present  Navier-Stokes 
solutions;  b)  other  results  exist  for  this  model 
and  hence  can  be  used  for  comparison;  and  c)  it  is 
a relatively  uncomplicated  model  to  compare  results 
for  hot  and  cold  flows. 

For  the  case  of  unsteady  flow,  basically  the 
same  flow  conf iquration  was  used  except  that  the 
flow  conditions  at  the  inlet  now  were  varied  in 
time  instead  of  keeping  them  constant.  In  the 
case  of  oscillations  in  the  .y  component  of  velocity 
at  the  inlet,  the  flow  is  off  course  not  a parallel 
mixinq  flow  although  the  time  averaqed  value  of  v 
velocity  was  zero.  Both  the  sinusoidal  and  sine 
square  fluctuations  at  the  inlet,  firstly  in  the 
F^  stream  and  then  in  both  the  streams  are  applied. 

III.  Analysis 

A.  Equations 

Each  of  the  nine  HF  vibrational  levels  is 
treated  as  a separate  species  interlinked  by  chem- 
ical and  vibrational  relaxation  reaction.  In 
addition  to  these  nine  HF  vibrational-level  species 
the  reactants  H,  H,,,  F,  F,,  and  an  inert  gas  (dil- 
uent) He,  make  a total  of  14  different  species 
leading  to  14  species  continuity  equations.  These 
equations,  in  addition  tothe  Navier  Stokes  equa- 
tions viz.  a)  qlobal  continuity  equation,  I ) x 
momentum  equation,  c)  > momentum  ec.iation  and  •') 
energy  equation  yield  a total  of  18  equations. 
Non-dimensional ized,  and  in  conservation  form, 
they  are 
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where  h is  the  characteristic  width  of  the  nozzle 
combination  as  shown  in  Fiq.  1,  subscript  r refers 
to  a characteristic  reference  value,  and  the  primes 
refer  to  dimensional  values. 
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Equation  ( 3 ) - ( 7 ) can  also  be  written  as 
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where  U,  F,  and  G are  vectors  with  18  components 
each,  containinq  the  terms  in  the  curly  brackets 
(6)  in  Eis.  (3)-(7).  Vector  R contains  the  riqht-hand 
sides  of  Eqs.  ( 3 ) - ( 7 ) . 


The  remaininq  unknowns  in  Eqs.  ( 3 ) - ( 7 ) are 
assumed  to  be  related  to  p,  T,  and  rk  in  the 
following  manner. 

Species  Internal  Enerqy  ek 
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where  ( orn, ) ^ ’s  the  heat  of  formation  of  species 


Calculation  of  Transport  Properties 

The  transport  properties  viz.  viscosity, 
thermal  conductivity,  and  diffusion  coefficients 
are  first  calculated  for  the  individual  species 
at  qiven  pressure  and  temperature  and  then  obtained 
for  the  complete  mixture  for  qiven  concentrations. 

Mixture  viscosity:  Viscosity  for  each  species  can 
be  qiven  by* ■ 
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Here  ' ' is  a function  of  , / kT  and  is  related 


The  previous  equations  are  nondimens ional ized  as 
fol lows 
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where  T*  = E/kT. 

The  vicosity  of  the  mixture  u,  is  then  obtain- 
ed using  the  Wilke  estimation  method  for  gases  at 
low  pressure^. 

6 6 y. 

u = I lnj/[l  ♦ E <t>i i („  )1  (13) 
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where 

c12/k=[(El/k)  x (e2/k)],/2 

variation  of  Tip  with  kT/t^,  is  available  in  a tab- 
ular form. 


where 
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The  196  (14  x 14)  binary  diffusion  coeffi- 
cients were  used  to  yield  multicomponent  diffusion 
coeff ic ;ents22  0 

im 


Thermal  conductivity:  The  thermal  conductivity  of 
an  individual  species  is  related  to  the  viscosity 
of  the  species  in  the  following  manner: 

21 

Monoatomic  species 
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= (15/4)(R/M)n  (14) 

21 

Diatomic  species 

A = *’(1 +0.88[(2/5)(Cp/R)-l])  (15) 

Thermal  conductivity  for  a mixture  of  gases 
at  low  pressure  is  given  by“ 

6 6 M 1/2  y. 

k * ^ (M f] 

jl»i 

6 6 y. 
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where 

V * A1(l/[H0.35(Cpi/R)-2)]) 

for  H^.Fg  and  HF  (diatomic  species) 

a id\*  * a.  for  H,F,  and  He  (monoatomic  spe- 
iei),  1 

also  x.**  = A ^ -xf  for  H2,  F?,  and  HF, 

* o for  H,  F,  and  He, 

"ij  * <W2)- 


Diffusion  coefficients:  Treating  each  vibrational 
level  of  HF  as  a separate  species,  the  binary  dif- 
fusion coefficients  are  first  obtained  for  all  14 
species;  these  binary  diffusion  coefficients  are 
independent  of  their  concentrations. 22 

Di j » O.OOieSST^^fM^MjJ/M.M^’^/po.^ap  (17) 

where  p is  the  pressure  in  atmosphere,  T is  tem- 
perature in  "K 

O 

j * (oj  * Oj )/2  in  A 
is  the  thermal  collision  Integral 


Species  diffusion  velocities  (or  diffusion 
mass  flux)  is  related  to  concentration  gradient 
by  Fick's  law. 


Chemical  Production  Term  1- : 

The  species  production  term  was  calculated 
explicitly.  Treating  each  vibrational  level  of 
HF  as  a separate  species,  100  elementary  reactions 
were  obtained  from  the  reactions  given  in  Table 
Vlll  Ref.  23.  These  reactions  involve  chemical 
pumping,  dissociation,  and  V-V  and  V-T  vibrational 
transfers.  The  vibrational  levels  of  H2  were  not 
treated  separately,  unlike  Ref.  23. 

The  equilibrium  constants  as  a function  of24 
temperature  were  obtained  from  the  JANAF  Tables ^ 
for  the  dissociation  reactions.  For  the  chemical 
pumping  and  V-V  and  V-T  transfers,  they  were 
assumed  to  be  of  the  form 

Keq  - exp[-(Ev-Eo)/RT] 

where  the  reaction  is  HFfOj^HFfv)  and  Ey  - Eq  is 

the  energy  of  the  vth  level  above  the  around  state. 
The  backward  rate  constants,  kp  were  then  obtained 

from  k.  = K /k,,  where  k,  is  the  forward  rate  con- 
d ea  t t 

stanr  obtained  from  Ref.  23. 

Therefore,  the  finite-rate  chemically  react- 
ing model  consists  of  100  elementary  reactions 
involving  14  species.  Hence,  14-species  rate 
equations  are  written,  one  of  which  is  given  below 
for  an  example. 

£ [HF(l)]=[k.31[H][F]-k3,[HF(l)])[M6] 

+(K4a[F][H2]-3k4a[HF(l)][H]) 

+(K5b[H][F2]-k_5b[HF(l)[F]) 

+(k_6a[HF(0)]-k6al[HF(l)])[M7] 


4 


/i>  5 mass  fraction 


♦(k6a2[HF(2)]-k6a2[HF(l)])[M7] 

+(k-6bltMF(°)]-k6bllHF(1 )])tM8] 
+(k6b2[HF(])]-K.6b2[HF(l)])[M8] 

*(k-6cItHF(0)3_k6bttHF(1 )])tM9] 
♦{k6c2tHF(2)]-k6c2[HF(l)])tM9] 

+(k-6cnCHF(0)]-k6ditMF(,)])[Mio] 

♦(k6d2[Hr(2)]-k6d2[HF{l)]KM10] 

♦ (k  6f1tHF(0)]-k.6f,tHF(l)])[H5] 


Cv  * lHF(v ) 


Q(v)  * E (2J+1  )exp-[E(v,J)1ic/kT]  (23) 

j-0 


* rotational  partition  function 
E(v,J)  = 20.9SJ(J+1)  - 0. 796J( J+l ) ( v+l/2)cm"'  (24) 

! vibrational -rotational  energy 


The  averaqe  integrated  gain  in  the  y direction 
can  be  obtained  as 


G 


v . J 


1 

h 


h 

J 

o 


It  . 

V*J 


dy 


(2b) 


♦(k6f2[HF(2)]-k_6f.2[HF(l)])[M5] 


C.  Boundary  Conditions 


*(‘.6ql[HF(0)]-k6g,[HF(l)])[M4) 

*(k6q2[HF(])]-k_6g2[HF(l)])[M4J 

♦2(k-7a,[HF(0)][HF(2)]-k7a[HF(l)]2) 

+(k7a[HF(2)]2-<_7d2[HF(l)][HF(l)]2) 
+(k_7b[HF(3)][HF(0)]-k7bfHF(l)][HF(  )]) 
+(k7b[HF(2)][HF(3)]-k_7b2[HF(l)][HF(4)]) 
Mk_7cl[HF(0)](HF(4)]-k7c[HF(l)][HF(3)]) 
♦(k7c[HF(2)][HF(4)]-k.7c2[HF(l)][HF(5)]) 


ft  comment  is  made  on  the  boundary  conditions 
at  y = 0 and  y = h.  Symmetry  conditions  hold  at 
these  boundaries,  as  can  be  seen  from  Fig.  1,  which 
mode  Is  a segment  of  the  multinozzle  flow  character- 
istic of  many  chemical  lasers,  i.e. 


•u  .-  2P  = 51  = M , v = o 
ay  ay  ay  ay 

at  y = 0 and  h.  In  the  present  f ini te-di fference 
scheme,  the  reflection  principle  is  used;  this  is 
an  accurate  representation  of  boundary  conditions 
on  a line  of  symmetry,  i.e. 


VI  *3-1 


■Vi 


V...  = -V.  , etc. 
J+l  J-l 


where  j is  an  index  in  y direction  and  j lies  on 
the  boundary  itself. 


Mk-7d,(HF(0(]tHF(5)]-k7d[HF(l)][HF(4)]) 
♦(k7a[HF(2)[HF(5)-k_7d  [HF ( 1 )][HF (6)]) 

B.  Gain  Calculations 

For  the  low-pressure  levels  that  exist  in  the 
mixing  region,  Doppler  line  broadening  is  assumed. 
Rotational  equilibrium  is  assumed  to  exist,  and 
only  the  P-branch  laser  transitions  are  considered 
Then  the  optical  small  signal  gain  coefficient  at 
line  center  has  the  value  of'^ 


VJ  = A(Cv+l  * vCv> 

(20) 

where 

a * 

(21) 

3(2kHHf./Kc)T/2Q(v+l)T'22 

exp  - 

\ , Otl+Ji  exp  - [E(v,J)-E(vH.J-l)hc/kT  (22) 
Q(v) 


IV . Numerical  Solution 

A time  dependent  technique  patterned  after  the 
MacCormack’5  approach  was  employed  to  qenerate 
the  steady-state  solution. 

Equation  (8)  can  be  written  in  the  following 

form. 


If  the  distribution  of  flowfield  variables  is 
specified  at  any  instant  of  time  (say  n),  the 
spatial  gradients  of  aF/ax  and  aG/ay  can  be  calcu- 
lated at  every  grid  point  by  finite  differences, 
and  the  vector  K can  be  explicitly  computed  from 
the  known  temperature,  rate  constants,  and  concen- 
trations. In  turn,  the  time  derivative  aO/at  can 
be  computed  from  Eq.  (26).  This  allows  the  compu- 
tation (in  principle)  of  new  values  at  time  (n+1) 
by 

Un+1  * 0"  + (|5)  At  (27) 

n 

The  components  of  the  new  vector  Gn+'  then 
specify  all  the  properties  at  a grid  point  one 
step  ahead  in  time.  Steady  state  is  essentially 
attained  when  aO/at.  approaches  zero.  However, 


Eq.  (27)  is  of  first-order  accuracy  only.  In  con- 
trast. NacCormack  uses  a predictor-corrector  method 
of  second-order  accuracy.  It  involves  the  genera- 
tion of  intermediate  predicted  values  at  time 
(n*l)  via  Eq.  (28).  Then  these  predicted  values 
are  used  again  in  the  conservation  equation  in  a 
"corrector"  fashion  to  ob’ain  values  at  time  (n+2). 
Averaging  of  these  two  steps  leads  to  a higher 
accuracy  at  time  (n*l)  as  shown  in  Eg.  (30).  This 
process  as  applied  to  mixing  flows  is  described  in 
more  detail  in  Ref.  26. 


° ♦ Rl  ,\t 

:iy  Jn 


* Hui  At 


Here  n refers  to  the  time  step  and  the  title  refers 
to  intermediate  values. 

The  vectors  T and  G contain  spatial  gradients 
of  temperature  and  velocities  directly  and  of  spe- 
cies densities  indirectly  through  the  diffusion 
velocities.  In  turn,  spatial  gradients  of  vectors 
P and  G themselves  need  to  be  taken  to  solve  Egs. 
(28)  and  (29).  Thus,  effectively,  spatial  gra- 
dients of  the  thermodynamic  guantities  need  to  be 
taken  to  a larqest  order  of  two.  This  offers  var- 
ious combinations  of  forward,  backward,  and  central 
differencing  schemes  to  compute  the  spatial  deriva- 
tives, a detailed  discussion  of  which  appears  in 
Ref.  14.  In  the  present  paper,  the  following 
differencing  scheme  was  used;  it  is  shown  in  Ref. 

14  to  be  the  most  appropriate  scheme  for  the  mixing 
flows  of  present  interest. 

Forward  differencing  is  used  on  the  P and  G 
vectors  to  obtain  the  intermediate  properties 
(predictor  step),  e.g. 


. r(j.j,k)  - P(2.j.k) 


Backward  differencing  is  used  for  the  corrector 
step.  e.g. 

(i  2 k)  * G^i,l,kJ  / j 

,iy  ' • \y  ' 


However,  a reverse  combination,  viz.  baikward 
differencing  for  predictor  step  and  forward  dif- 
ferencing for  the  corrector  step,  is  used  for  the 
gradients  of  T,  u,  v,  o,  and  e.  . This  is  termed 

the  "Modified  MacCorma  k"  approach  in  Ref.  14. 

The  time  step  used  to  advance  the  solution  in 
time  is  the  ijjinimum  of  the  time  steps  due  to  the 
CFL  criteria-'  and  the  characteristic  relaxation 
time.  Here,  the  CFL  criteria  are  given  by  the 
minimum  of 

At  * Ax/(u*a)  and  At  - Ay/(y*a) 

28 

whereas  the  chemical  relaxation  time  is  given  as 


A 9 x 9 grid  size  (81  mesh  points)  is  used. 


V.  Results  and  Discussion 
A.  Steady  Flow 

Solutions  are  obtained  for  the  flow  depicted 
in  Fiq.  1.  At  time  t * 0,  the  initial  conditions 
are  rather  aribtrarily  chosen  as  constant  proper- 
ties in  each  of  streams  1 and  2 (same  as  the  up- 
stream boundary  conditions  at  the  nozzle  exits 
held  fixed  with  time). 

Consider  the  point  defined  by  x/h  10.0  and 
y/h  - 0.375.  The  time  history  of  the  static  tem- 
perature at  this  point  is  given  in  Fig.  2.  Cold- 
-flow  calculations  (the  fluid  dynamics  without 
the  chemical  reactions)  are  made  through  a nondi- 
mensional  time  of  24  to  allow  the  two  streams  to 
partially  mix.  After  this  time,  hot-flow  calcula- 
tions (the  fluid  dynamics  fully  coupled  with  the 
chemical  kinetics)  are  made.  Note  from  Fig.  2 that 
the  switch  to  a hot  flow  causes  a discontinuous 
increase  in  dT/dt,  and  a subsequent  approach  to  a 
hiqher  steady-state  temperature.  For  comparison, 
the  purely  cold-flow  case  is  carried  out  to  a 
steady  state,  as  also  shown  in  Fiq.  2.  Note  that 
the  combined  effect  of  "domical  reactions  and 
vibrational  relaxation  .eao  to  a 35T  increase  in 
static  temperature  in  comparison  to  the  purely 
aritifical  cold-flow  case. 

Steady-state  profiles  of  velocity,  pressure, 
and  temperature  with  respect  to  y/h  are  given  in 
Fiqs.  3.  4,  and  5,  respectively.  In  each  figure, 
results  for  the  longitudinal  stations  x/h  0,5 
and  10  are  given.  Also  in  each  figure,  the  solid 
and  dashed  lines  correspond  to  hot  and  cold  flow, 
respectively.  Note  from  Fig.  3 that,  in  contrast 
to  temperature,  the  velocity  profiles  downstream 
of  the  nozzle  exits  show  little  influence  due  to 
chemical  reactions  --  a result  that  is  almost 
classical  in  most  high- temperature  gasdynamic  pro- 
blems. Also  note  that  the  slower-moving  stream  of 
H.,  is  accelerated  more  than  the  faster-movino 

C 

stream  of  F and  F,  is  retarded,  due  to  the  mixing 

process.  The  pressure  variations,  which  are  plot- 
ted on  an  expanded  scale  in  Fig.  4,  show  virtually 
a constant  pressure  in  the  y direction,  except 
for  a small  variation  in  the  middle  of  the  mixing 
zone.  For  hot  flow,  there  is  an  adverse  pressure 
gradient  in  the  flow  direction;  in  contrast  the 
net  effect  of  the  purely  fluid  dynamic  mixing 
appears  to  be  a slight  favorable  pressure  gradient, 
at  least  for  the  first  10  nozzle  heights  down- 
stream. The  corresponding  temperature  profiles 
are  shown  in  Fig.  5.  Comparing  the  cold-flow  pro- 
files. it  is  apparent  that  the  viscous  shear  action 
causes  a larqer  local  temperature  rise  in  the  mix- 
ing region  at  x/h  * 5.0  than  at  x/h  * 10  because  y 
gradients  of  velocity  are  larger  near  the  nozzle 
exits.  However,  in  going  from  x/h  - 5.0  to  x/h 
10.0,  the  mean  temperature  rise  is  more  over  the 
whole  cross  section  whereas  the  peak  temperature 
of  the  faster  streams  seems  to  drop  because  the 
interface  streamline  would  bend  toward  the  slower 
stream.  The  average  pressure  also  si ight  ly decreases 


* 


as  the  flow  moves  downstream,  as  noted  in  Fig.  4. 

In  contrast,  the  hot  flow  profiles  clearly  show  an 
almost  discontinuous  increase  in  the  temperature, 
a fact  already  noted  from  Fig.  2. 

figures  6 and  7 illustrate  density  profiles 
of  various  HF  vibrational  levels  at  x/h  5 and  10, 
respectively.  The  growth  of  the  reaction  zone  can 
he  seen  clearly.  Total  population  inversions  exist 
between  vibrational  level  increases  as  the  flow 
moves  downstream,  the  small  signal  gain  does  not 
necessarily  increase  since  the  absolute  difference 
between  the  densities  of  HF(0),  HF ( 1 ) , and  HF(2) 
does  not  always  increase. 

The  calculated  growth  of  the  reaction  zone,  as 
defined  by  the  region  where  is  greater  than 

10  of  the  maximum  value,  is  shown  in  Fig.  8.  This 
figure  is  taken  from  Of.  12,  and  contains  experi- 
mental data  taken  at  TRW.  The  present  results  are 
marked  on  Fig.  8 and  show  the  same  laminar  varia- 
tion as  the  experimental  data. 

The  existence  of  population  inversions  and 
hence  laser  action  is  best  seen  in  Fig.  9.  where 
densities  of  various  HF  vibrational  levels  at 
y/h  - 0.375  are  plotted  with  respect  to  x. 
Inversions  exist  between  the  0-1  levels  and  1-2 
levels,  which  yields  the  small-signal  gain  .«  , 

as  calculated  by  Eg.  (20).  These  gains  are  shown 
in  Figs.  10  and  11.  Figure  10  shows  the  variation 
of  the  small  signal  gain  with  respect  to  x at 
y/h  » 0.375  for  the  vibrational  level  transition 
1-0  for  various  rotational  levels  (only  P-branch 
transitions  are  assumed).  Figure  11  shows  similar 
results  for  the  2-1  transition.  The  values  of 
gain  and  the  spatial  extent  of  the  lasing  region 
as  indicated  on  Fig.  10  and  11  are  typical  of 
conventional  HF  chemical  lasers,  as  obtained  from 
Refs.  18  and  19. 

A few  words  about  stability  and  convergence 
are  in  order.  The  present  time-dependent  calcula- 
tions smoothly  and  regularly  approach  a steadv- 
-state  solution  as  long  as  the  requisite  stability 
criteria  are  followed;  i.e.,  the  solutions  are 
stable  as  long  as  the  time  increment  is  less  than 
the  CFl  and  chemical  relaxation  times.  With  regard 
to  convergence,  the  question  can  be  asked:  Is 
enough  accuracy  obtained  with  the  present  9x9 
grid,  which  at  first  glance  appears  rather  coarse? 
An  answer  is  given  in  Fig.  12.  Here,  the  final 
steady-state  temperature  at  x/h  * 5 and  y/h  = 0.75 
is  given  for  three  different  grid  sizes:  5x5, 
9x9  and  13  x 13.  It  appears  that  a 9 x 9 grid 
is  sufficiently  accurate,  and  that  a further 
definition  by  more  grid  points  is  unnecessary. 

This  is  totally  consistent  with  time- dependent 
solutions  of  of  other  problems  1 .29-31  where  suffi- 
cient accuracy  has  been  obtained  with  seemingly 
very  coarse  grids.  Apparently,  the  time-dependent 
mechanism  is  "self-correcting"  at  each  time  step, 
allowing  the  physics  contained  in  the  conservation 
equations  to  bear  more  strongly  and  accurately 
at  each  grid  point.  The  philosophical  point  not 
withstanding,  experience  has  clearly  proven  that 
time-dependent  solutions  require  fewer  grid  points 
than  might  be  expected  for  steady-state  analyses. 
The  present  results  are  a case  in  point.  As  long 
as  the  cell  Reynolds  number  in  the  y direction  is 
on  the  order  of  unity,  the  gradients  are  adequately 
accounted  for,  as  in  the  present  calculat ions . 


lhe  present  Navier-Stokes  calculations  are 
compared  with  the  results  of  King  and  Mirels'9  and 
1AMp'h  in  Figs.  13  - 15.  Both  Refs.  18  and  19  are 
boundary-layer  solutions.  King  and  Mirels'^  assume 
two  semi-infinite  streams;  hence,  the  pressure  is 
assumed  constant  in  both  the  x and  y directions. 

The  LAMP  program  calculates  a pressure  gradient  in 
the  x direction  by  means  of  a quasi-one-dimensional 
heat-addi t ion  approximation,  but  assumes  constant 
pressure  in  the  y direction.  The  present  Navier- 
-Stokes  calculations  allow  pressure  gradients  in 
any  direction.  The  authors  express  their  extreme 
appreciation  to  Dr.  Walter  dlowacki  of  the  Naval 
Surface  Weapons  Center  (White  Oak,  Md.)  for  his 
running  of  the  LAMP  code  for  comparison  with  the 
present  work. 

The  caseshownin  Fiqs.  13-15  is  a strongly 
reacting  mixture,  with  conditions  at  the  nozzle 
exit  of  (see  Fig.  1 for  the  nomenclature)  p - 5 
Torr,  T.  H0°K,  T,  400°K,  u.  1400  m/sec, 

I t I a 

u.,  2140  m/sec.  In  stream  1,  . 1.47  x 10 

i H? 

kg/m'.  In  stream  2,  , F 1.1  x 10  kg/m  and 
-4  1* 

, 5.7  x 10  kg/m  . This  mixture  is  so 

e 

strongly  reacting  that  a large  adverse  pressure 
gradient  is  produced  in  the  flow  firection;  note 
from  Fig.  13  that  both  the  present  calculations  and 
LAMP  predict  approximately  a factor  of  five  pres- 
sure increase  in  a distance  on  the  order  of  a cen- 
timeter. In  fact,  the  changes  are  so  severe  that 
the  LAMP  calculation  experienced  some  spurious 
wiggles,  and  then  blew  up  beyond  0.8  cm.  (However, 
this  in  no  way  reflects  on  the  viability  of  LAMP, 
because  no  subsequent  efforts  wpre  made  to  adjust 
grid  size,  etc.,  to  attempt  to  overcome  this  beha- 
vior.) Nevertheless,  the  comparison  be’ween  LAMP 
and  the  present  calculations  is  reasonable;  in 
contrast,  the  constant  pressure  assumption  of 
King  and  MirelslR  is  not  valid  for  this  case. 
Similar  comparisons  for  T and  HF  variations  are 
given  in  Figs.  14  and  15,  respectively. 

Return  to  the  conditions  given  in  Fiq.  1. 

Here,  the  densities  of  F,,  F and  H are  relatively 
small  --  a dilute  case.  'The  results  of  Figs.  2-11 
apply  to  this  case.  In  turn,  these  results  are 
compared  with  a calculation  from  LAMF  in  Fig.  16. 
The  profiles  of  HF ( 0 ) and  T are  shown  as  a function 
of  y at  x/h  » 10.  Note  that  the  LAMP  results 
predict  peak  temperatures  and  HF  densities  that 
are  about  20  and  70'  higher,  respectively,  than 
the  present  calculations . This  is  considered  to 
be  reasonalbe  agreement  in  light  of  the  different 
kinetic  rates  and  transport  properties  in  the  two 
programs.  The  present  results  calculate  detailed 
and  variable  transport  properties  at  each  point  in 
the  flow;  lamp  assumes  constant  Prandtl  and  Lewis 
numbers.  Also,  keep  in  mind  that  the  pressure 
gradients  in  the  two-dimensional  flow  appear  quite 
naturally  and  exactly  in  the  present  Navier-Stokes 
formulation,  whereas  l AMP  has  an  approximate  cal- 
culation of  pressure  gradient  and  only  in  the  x 
direction. 

A word  about  computer  time  is  also  in  order. 
For  the  present  Navier-Stokes  results,  computer 
times  of  30  and  90  min  on  a 0N1VAC  1108  for  the 
dilute  and  stronqly  reacting  cases,  respectively, 
were  required.  The  cost  averaqed  $350  per  run. 

This  is  about  three  times  more  expensive  than  LAMP 
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for  the  same  cases.  However,  keep  in  mind  that  the 
present  time-dependent  calculations  start  with 
initial  conditions  that  are  constant  but  different 
properties  in  streams  1 and  2,  i.e.,  no  mixing  of 
the  two  streams.  If  more  realist tc  initial  condi- 
tions are  fed  in,  the  computer  time  can  be  drasti- 
cally reduced.  In  particular,  if  several  different 
parametric  runs  are  stacked  together,  and  the  ini- 
tial conditions  for  one  are  taken  as  the  solution 
from  another,  the  running  times  should  be  compar- 
able to  or  better  than  LAMP.  This  behavior  is 
already  observed  in  time -dependent  calculations 
of  gasdynamic  laser  performance,^  where  by  stack- 
ing several  runs,  the  running  times  of  the  second 
and  remaining  run s are  one-fifth  that  of  the  first 
run.  In  this  fashion,  Navier-Stokes  calculations 
can  be  made  cost-effective;  therefore,  they  are 
not  inordinately  long  as  may  first  be  expected. 

B.  Unsteady  flow 

The  main  thrust  of  the  present  paper  is  to 
answer  what  happens  if  one  (or  more)  of  the  vari- 
ables at  the  inlet  in  one  of  the  streams  (or  both 
streams)  is  fluctuated  in  line.  Does  it  enhance 
or  inhibit  mixing  and  how  does  it  affect  gain? 

It  was  observed  that  a time-wise  sinusoidal 
variation  in  inlet  pressure  in  the  F,  stream  (and 

hence  density,  keeping  the  temperature  constant) 
about  a certain  mean  value  failed  to  improve 
mixing.  Correspondingly  no  better  gain  was  obser- 
ved. Variation  in  frequency  and/or  amplitude  of 
the  pressure  fluctuation  did  not  improve  this 
situation. 

Next,  in  order  to  vary  the  momentum  of  the 
incoming  flow,  the  x-component  of  velocity  was 
fl uctuated  about  a mean  val ue  of  2000  m/sec.  A 
noticeable  effect  of  mixing  was  observed  only  for 
angular  velcoity  about  3 x 10^  rad/sec  and  for  a 
rather  large  amplitude  of  +500  m/sec.  A plot  of 
velocity  vectors  in  the  flow  at  a particular  time 
indicated  the  presence  of  a large  induced  y-compo- 
nent  of  velocity  in  the  flow  due  to  coupling  with 
the  fluid  dynamics,  and  this  is  perhaps  responsible 
for  increased  mixing. 

Therefore,  all  subsequent  cases  were  run 
wherein  only  the  sinusoidal  variation  of  the  y- 
-component  of  velocity  was  intentionally  forced  at 
the  inlet.  This  yielded  a considerably  improved 
mixing  and  also  higher  values  for  small  signal 
and  in  particular  the  integrated  gain.  Many 
different  cases  were  simulated  of  which  a case 
with  frequency  ...  * „ x 10^  rad/sec  and  an  ampli- 
tude of  v = ♦ 500  m/sec  yielded  the  most  favorable 
results.  For  example.  Figure  17  illustrates  the 
transient  variation  of  the  F.,  density  at  a given 

point  in  the  flow  for  three  different  amplitudes 
in  the  fluctuation  of  the  inlet  y-component  of 
velocity  (v),  namely  \v  = + 5'-,  + 10/  and  + 25 
of  the  x-component  of  velocity  (u).  The  horizon- 
ta)  dashed  lines  are  the  F^  densities  obtained 
from  a previous  steady-flow  solution  for  the  same 
inlet  properties.  Note  that  for  \v  * + 25’  the 
unsteady  flow  results  in  a time-mean  F ? density 
that  is  above  the  steady-flow  value,  these  un- 
steady results,  when  integrated  over  time,  lead 
to  the  results  in  Figure  18.  Here,  the  flow-wise 
variation  of  time  averaged  density  of  Fp  is  shown 
for  various  cases.  This  shows  clearly  the  improve- 


ment in  mixing  due  to  the  fluctuating  of  v at  the 
inlet.  Gain  results  also  show  the  consequence  of 
the  enhancement  of  mixing.  A rather  remarkable 
factor  of  two  increase  in  the  integrated  gain  was 
obtained  for  the  case  outlined  above  (as  discussed 
later) . 

It  was  further  observed  that  oscillations  of 
v in  both  the  F^  and  H^  streams  at  various  differ- 
ent phase  angles  yielded  results  generally  similar 
to  or  inferior  to  the  fluctuations  in  only  one 
stream. 

It  should  he  emphasized  here  that  the  frequen- 
cy of  oscillations  was  selected  to  be  of  such  a 
magnitude  so  as  to  have  wavelength  of  the  order  of 
the  cavity  length.  Further  parametric  studies 
confirmed  this  hypothesis.  The  Strouhal  number  is 
of  the  order  recommended  in  Ref.  33. 

Figures  19  and  20  clearly  indicate  the  remark- 
able superiority  of  unsteady  mixing.  All  the 
results  here  have  again  been  time  averaged  over  a 
periodic  time  of  one  oscillation.  Figure  19  shows, 
clearly,  a larger  and  more  spread  out  production 
of  HF  for  the  case  of  v velocity  fluctuation  *han 
the  other  two  plotted.  Ability  of  the  oscillation 
to  populate  the  region  of  the  flow  for  y/h  0.5, 
which  did  not  occur  for  the  other  two  cases,  in 
general  produced  integrated  gain  larger  by  a factor 
of  two,  thus  making  it  possible  to  increase  gain 
by  a substantial  amount.  Figure  20  again  shows 
various  HF  levels  versus  x for  a constant  value  of 
y.  Here  again  a faster  mixing  and  faster  pumping 
is  clearly  seen  e.g.  at  x/h  - 2.5. 

2 

Sine  square  [sin  <•]  variation  in  the  y-compo- 
nent of  velocity  was  also  simulated,  thus  yielding 
a definite  value  for  the  time  averaged  mean  of  a 
fluctuation.  This  suggested  another  case  of  non- 
-parallel  mixing  of  two  steady  streams  at  an  angle 
(viz.  At'  - 14°)  w.r.t.  each  other.  Both  these 
cases  produced  inferior  mixing  gain  results.  The 
shock  waves  produced  in  the  later  case  accounted 
for  a rise  ir  temperature  in  the  lasing  region 
and  hence  taster  V-V,  V-T  deactivation  reactions, 
thus  impeding  gain. 

The  uniform  profiles  of  velocity  and  concen- 
trations of  F and  F,,  used  at  the  inlet,  although 

sufficient  and  proper  for  parametric  studies  at 
earlier  time,  need  to  be  modified  to  account  for 
no  slip  condition  and  catalytic  recombination  at 
the  nozzle  wall.  Also  a recirculation  region  is 
produced  near  the  inlet  due  to  finite  thickness 
of  nozzle  walls.  These  questions  are  being 
addressed  to,  presently,  using  the  Navier-Stokes 
equations  to  solve  the  flowfields.  In  particular 
the  Boundary  Layer  profile  in  the  x-component  of 
velocity  at  the  inlet  in  both  streams  has  already 
been  simulated.  The  final  integrated  gain  results 
are  shown  in  Figure  21  for  P-transiion  with  V,J 
1,1  for  various  different  aforementioned  cases 
including  the  boundary  layer  case.  The  decrease 
in  qain  due  to  the  presence  of  real  effects  like 
B.  L.  is  easily  seen.  Integrated  gain  is  plotted 
in  I iqs.  22  and  23-  In  Fig.  22  integrated  gain 
plotted  versus  non-dimensional  x-wise  location 
clearly  illustrates  the  increase  in  gain  exper- 
ienced due  to  the  unsteadiness  of  the  flow.  For 
the  case  of  transitions  from  vibrational  levels 
I to  0.  Again  Fig.  23,  which  is  a similar  plot 
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for  the  vibrational  level  transition  2 to  1 reem- 
phasizes the  same  point.  In  both  these  cases  the 
gain  is  calculated  from  the  various  HF  level 
species  densities  increased  over  a time  period  of 
the  osci 1 lation. 

As  a final  conment  on  the  fluid  dynamic  aspects 
of  unsteady  mixing,  note  that  there  is  a physical 
improvement  in  mixing  due  to  the  forced  interming- 
ling of  the  fluid  elements.  This  can  be  seen  from 
Fig.  24a  and  b,  which  shows  the  time  history  of 
the  particles  inserted  at  two  different  times  in 
the  flow.  The  nonuniform  particle  paths  clearly 
show  the  extent  of  material  interminql ing.  Those 
path  lines  were  obtained  as  part  of  the  computer 
results  by  means  of  a special  subroutine. 


V.  Conclusions 

The  present  paper  introduces  a "third  genera- 
tion" of  chemical-laser  analysis,  i.e,  Navier- 
-Stokes  solutions  for  the  flowfield  coupled  with 
the  detailed  chemical  kinetics  for  both  the  hot 
and  cold  reactions  of  HF  for  both  steady  and  un- 
steady flows.  The  results  for  the  steady  flow  are 
also  described  in  Ref.  35.  In  particular,  the 
present  results  show  the  following: 

1.  Navier-Stokes  solutions  for  supersonic 
diffusion  chemical  lasers  are  indeed  feasible; 
however,  computer  tines  equivalent  to  about  30 
min  or  longer  on  a UNIVAC  1108  are  required  for  a 
single  case.  By  stacking  cases  back-to-back  such 
that  the  initial  conditions  for  one  are  obtained 
from  the  solution  of  another,  the  net  time  per  case 
can  be  substantial ly  reduced. 

2.  The  major  potential  for  such  Navier-Stokes 
-olutions  is  in  the  analysis  of  recirculation  and 
separated  flow  effects  on  laser  performance. 

3.  In  a comparison  between  hot  and  cold  flows, 
the  chemical  reactions  markedly  affect  the  tempera- 
ture distributions,  but  have  little  effect  on  the 
velocity  distributions.  The  pressure  increases  in 
the  flow  direction  due  to  chemical  reactions,  an 
effect  to  be  expected  from  simple  analogy  with 
constant  area  heat  addition  in  supersonic  flows. 

In  contrast,  for  cold  flows,  the  longitudinal 
pressure  variation  is  reasonably  constant  and  may 
even  decrease  slightly. 

4.  The  growth  of  the  laminar  reaction  zone 
predicted  in  the  present  paper  compares  favorably 
with  experiment. 

5.  Navier-Stokes  calculations  have  the  dis- 
tinct advantage  that  the  two-dimensional  pressure 
gradients  appear  quite  naturally  and  exactly.  In 
strongly  reacting  cases,  the  proper  accounting  of 
these  pressure  gradients  are  absolutely  necessary; 
the  constant  pressure  boundary- layer  assumption  is 
not  adequate. 

6.  Considering  the  differences  between  the 
physical  properties  and  fluid  dynamic  modeling 
of  the  present  Navier-Stokes  analysis  and  LAMP, 
fairly  reasonable  agreement  is  obtained  between 
the  two.  However,  the  constant  pressure  results 
of  King  and  Mirels  differ  considerably,  and  are 
obviously  not  valid  for  cases  where  strong  pressure 
gradients  occur  in  the  real  problem. 


mm 


7.  It  is  the  first  numerical  solution  of  the 
full  Navier-Stokes  equations  with  fully  coupled 
kinetics  for  both  the  cold  and  hot  reactions  for 
HF  and  multicomponent  diffusion  for  the  case  of 
inherently  unsteady  supersonic  mixing  flows. 

8.  It  shows  that  the  integrated  gain  for  HF 
chemical  lasers  can  be  increased  substantial ly 

by  intentionally  fluctuating  the  y-component  of 
velocity  at  the  inlet.  This  leads  to  better  mixing 
and  also  higher  gain  in  general. 

Real  effects  like  presence  of  bounary 
layer  in  velocity  at  inlet  impeded  mixing  and 
hence  reduces  integral  gain. 

The  authors  recognize  that  the  density  varia- 
tions created  in  the  cavity  by  the  fluctuations 
will  have  an  adverse  effect  on  the  overall  quality 
of  the  laser  beam.  However  the  purpose  of  this 
study  has  been  to  examine  the  fluid  dynamic  and 
chemical  kinetic  consequences  of  such  fluctuations. 
The  effect  on  beam  quality  is  beyond  the  scope  of 
the  present  paper. 
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Geometric  model  and  given  cavity  inlet  conditions. 
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Figure  2 

Transient  variation  of  static  temperature  at  x/h 
1 .0  and  y/h  = 0.375. 
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Figure  3 

Steady-state  velocity  profiles  at  various  x-wise 
locations. 
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Figure  4 

Steady-state  pressure  profiles  at  various  x-wise 
locations. 
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Figure  5 

Steady-state  temperature  profiles  at  various  x-wise 
locations. 
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Steady-state  density  profiles  for  HF(v);  x/h  = 5.0. 
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Figure  7 

Steady-state  density  profiles  for  HF(v);  x/h  10.0. 
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Figure  8 

Growth  of  the  reaction  .rone;  comparison  of  present 
calculations  with  experiment  (see  Ref.  1?). 
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Variation  of  small-signal  gain  with  flow  direction 
2-1  transition,  y/h  0.375. 
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Figure  8 

Variation  of  HF(v'  uh  the  flow  direction;  y/h*0.375. 
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Temperature  vs.  grid  size;  a test  of  convergence 
propert ies . 
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Figure  13 

Pressure  vs  longitudinal  distance;  comparison 
between  present  results  and  methods  of  Refs.  18 
and  Id.  Srongly  reacting  case. 
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Temperature  vs.  longitudinal  distance;  comparison 
between  present  results  and  methods  of  Refs.  18 
and  19.  Strongly  reacting  case. 
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HF(0)  density  vs  longitudinal  distance;  comparison 
between  present  results  and  methods  of  Refs.  18  and 
19.  Strongly  reacting  case. 


Figure  lo 

HF(0)  and  temperature  profiles;  comparison  between 
present  results  and  methods  of  Ref.  IS.  Dilute 
case. 
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Figure  22 

Axia)  variation  of  integrated  gain;  1-0  transition 
steady  and  unsteady  flows. 


Figure  23 

Axial  variation  of  integral  gain;  2-1  transition. 
Steady  and  unsteady  flows. 
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PARTICLES  INSERTED  AT  PHASE  ANGLE  ♦■0  of»'0 


PATHLINES  OP  INERT  PARTICLES 
INSERTED  AT  INLET  AFTER 
REPEAT  IVE  STATE  IS  REACHED 


Figure  24 

Path  lines  of  given  fluid  elements  in  the  fluctua- 
ting flow.  Sinusoidal  fluctuations  in  v velocity 
at  F2  stream  inlet  for  inert  particles  inserted  at 

two  different  phase  angles,  w ="  x 105  rad. /sec. 
and  v = + 25%. 
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